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Chiral Lagrangians at finite density 
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The effective SU(2) chiral Lagrangian with external sources is given in the presence of non- 
vanishing nucleon densities by calculating the in-medium contributions of the chiral pion-nucleon 
Lagrangian. As a by product, a relativistic quantum field theory for Fermi many-particle systems at 
zero temperature is directly derived from relativistic quantum field theory with functional methods. 



PACS numbers: 12.39.Fe, 24.85.+p, 13.75.L, 13.75.Gx 

1. In the limit of massless up and down quarks the 
QCD Lagrangian is symmetric under the chiral group 
SU(2)l x SU(2)r. One assumes that this symme- 
try is spontaneously broken to the diagonal subgroup 
SU(2)l+r giving rise to the appearance of 3 massless 
Goldstone bosons which finally acquire small masses due 
to the non- vanishing mass of the u and d quarks. 

This symmetry breaking scenario constrains so much 
the interactions of the Goldstone bosons that the QCD 
Green functions can be calculated at low energies as an 
expansion in powers of momenta and quark masses. This 
is known as Chiral Perturbation Theory 

The extension of the theory to the case of low temper- 
ature at zero density was considered in ref. || . In this ar- 
ticle we study the case of small densities at zero temper- 
ature and derive the corresponding chiral Lagrangian by 
calculating the in-medium contributions due to the chiral 
pion-nucleon Lagrangian |^,^| with functional methods. 
Although we focus our treatment to QCD, the relativistic 
many-body formalism here deduced for Fermi systems 
can be applied to processes governed by other dynamical 
theories, as the traditional non-relativistic zero temper- 
ature many-body quantum theory which stresses 
the diagrammatic approach. Compared with standard 
quantum field theory at finite temperature T [^|J^] in the 
grand canonical ensemble, one avoids the use of unknown 
chemical potentials which themselves have to be calcu- 
lated in terms of the many-body forces. The former is 
accomplished by following quantum field theory at T = 
considering directly the change of the ground state from 
the vacuum to one with finite fermionic densities. In the 
same way one also avoids the non-trivial T — > limit due 
to the so called anomalous diagrams flfj| , p| . The price 
to pay is to rely on the adiabatic hypothesis in order to 
determine the interacting ground state from that of the 
free case by turning on the interactions adiabatically. 

2. Let us take first the case of symmetric and unpo- 
larized nuclear matter, the extension of the formalism 
to the asymmetric and polarized case is straightforward 
and will be shown below. In the following we take the 
Heisenberg picture and, following closely scattering the- 
ory p3[, we consider two ground states l^out) an d l^in) 
which, under the action of any time dependent opera- 
tor at asymptotic times t — > ±oo, respectively, behave as 
two symmetric Fermi seas of free protons and neutrons. 



The Fermi seas are filled up to the corresponding bary- 

JV 

onic density, JJa^(p n )|0), where the label n includes also 
the spin and isospin indices, N is the number of momen- 
tum states inside the Fermi sea with Fermi momentum 
kp = (37r 2 p/2) 1 / 3 , p is the total nuclear density and |0) 
is the vacuum. Our objective is to evaluate the gen- 
erating functional Z[v, a, s,p] in the presence of vector 
axial a p , scalar s and pseudoscalar p external fields 
by working out the transition amplitude (^outl^in)^i 
where the label J just indicates the presence of the afore- 
mentioned external sources. In this way by taking func- 
tional derivatives of Z[v, a, s, p] with respect to the exter- 
nal sources one evaluates the in-medium QCD connected 
Green functions (space-time averages at finite density of 
the quark currents coupled to the v, a, s and p sources). 
To do this we consider the effective chiral Lagrangians 
C = C^iT + Cm + AavW + ... with increasing number of 
pairs of nucleon fields ^>{x). We first restrict ourselves to 
the term with no nucleon fields and to that contain- 
ing two of them — ijj(x)D(x)'ip(x), together with the 
previous external fields. We will discuss later a way to 
include pcrturbatively the contributions of Lagrangians 
with higher number of nucleons by considering them to 
arise from bilinear vertices through the exchange of an 
arbitrary heavy particle. Indeed, although we are talking 
about CHPT, the only thing that matters for the follow- 
ing derivations is that is bilinear in the fermions. 
Consider now the transition amplitude for the ground 
states from t — > — oo to t — > +oo in the presence of 
the previous external sources together with Grassmann 
sources -q and rf, coupled to the nucleon fields: 

(fioutlfiinUwt = J[dU}[dvh][d^}(n out \i:(+^)) x 
e iSdx[c^+i>D^i>+^v] (^(_oo)m n ) (i) 

with the pion fields described by the 2x2 unitary matrix 
U. 

The ground state functional (^>(±oo)|fi out ) can be 

expressed in terms of that of the vacuum by writing the 
a(Pn)out operators as a function of ip{x) and its time 



derivative ip(x): 



a(Pn)out = lim 



E(p) 



t— >±oo 2m 



N 



e lE ^ 1 /dxe-*P x u(p„)[V^) 



1 



+ E^j^(x)}, ( 2 ) 

where as usual in scattering theory for t ±00 
the matrix elements are calculated as if there were 
no interactions. In the former expression, E(p) is 
the energy of the nucleon with three- momentum p, 
?7i7v is the nucleon mass and tt(p n ) is a Dirac spinor. 
Expressing tp(x) in terms of ip(x) (one way is us- 
ing the Dirac equation), taking into account that 

^ x ySdy^(y)^(y) = _j_^ e » f <W ( a )</>(y) ^ the anak) _ 

gous expression for t/^x), and substituting all that in 
eq.(0), one obtains for (fi utl fi in) J,vrt : 



lim [\dU](detD)e 

t — ► +00 / 



if dxC-n 



f N 

n 



E(Pn) 

2m,N 



u(Pn) 



d 



E( P n) 



(7° '^J + ^Vv) 



N 



-ifdxJdy^(x)D- 1 {x,y)j a r 1 {y) j~| E{q m ) 



5v(Vm) E{q m ) dyfa Sr)(y m ) 



2m n 

7 




7 fc 7° + i 



5v(ym) 



7°u(qm) 



(3) 



where the integration over the nucleon fields and con- 
jugate momenta is also done. Furthermore we define 
x n = (t,x„), y m = (t',y m ), p n = (E(p n ),p n ) and anal- 
ogously for q m . The action of the spatial derivatives can 
be readily taken into account by integrating by parts. In 
this way, they only act on the corresponding exponen- 
tials giving rise to three-momenta factors that can be 
further simplified by applying the Dirac equation on the 
Dirac spinors. In this way we have: 



"(Pn) 



h E{pn) E{pn) 



E( P r, 



tVp«) 



2m N 

2m,N 



HPr. 



Applying the previous results to eq.(0) it simplifies to: 



lim J [dU](det D)e^ dxC ^ (^Yljdx 



5rf(x n ) 

7 



,—if dxj dyrf (x)D 1 (x ,y)~f° ri(y) 



N 



u{q m ) 



After acting with the left derivatives 5 /5rf(x n ) on the 
exponential depending on the Grassmann sources, the 
former expression can be recast as: 

^ lim j [dU]{det D)e^ dxC ^ ( ]J f dx n e 1 "" 1 " f dz n , 

N 



This result is equal to e lZ ^ v ' a ' s ' p ^ when rj — > and 77+ — » 
O.Q Hence we can equalize to 1 the second exponential 
from the left and the remaining right derivatives and 
T](z n ) sources have to be paired in order to finish with a 
non- vanishing result. Thus we have: 



e iZ = t lim / [dU] (det D)e^ dxC ^ J2 e ^f[[ dx " 

Jd Xan e^u^p^D-^x^x^fe'^^uip^) (4) 

with e(cr) the signature of the permutation a over all the 
indices -momenta, spin and isospin. In order to continue 
let us write the operator D{x) = Dq(x) — A(x) with 
Dq{x) = ij^d^ — m n the Dirac operator for the free 
motion of the nucleons. On the other hand, the operator 
A(x) is completely general although in our case at hand, 
CHPT, it is subject to a chiral expansion of powers of 
soft three-momenta and quark masses. Furthermore, let 
us note that: 

lim / dxe ipx uUp)Do 1 (x,x') = -ie ipx 'u(p) 

i-^ + OO J 



lim dxDv"{x',x) 1 [> e- lvx u{p) = -ie- lpx u{p) (5) 

This result can be easily obtained by writing Dq 1 (x, x') 
in four-momentum space and then performing the in- 
tegral over the temporal component of the momentum 
taking care of the imposed limits. For instance, let us 
take the first of the previous equations. Then we have: 

lim [ dxe ipx uUp)Do 1 {x,x') = 

i— »-f oo J 

dR e- lR ( x - x '^(R+m N ) 



t— >+oo 



lim / dxe lp V(p) 



(2tt) 4 R 2 -m% + ie 



'Without taking the limit r), ry — > 0, the generating 
functional contains baryonic sources and taking differential 
derivatives with respect to them one could directly evaluate 
in-medium baryonic Green functions, e.g. nucleon propaga- 
tors. Nevertheless, since for our present purposes the baryons 
in the medium constitute just a background, we do not con- 
sider this case any further. 



2 



with e a positive infinitesimal. Exchanging the order of 
the integrations, the spatial one gives rise to (27r) 3 <5(R — 
p) which fixes R. As a result one has: 



lim 

t — >+oo 



— M t( p ) e 6 J. (BP 

2ir i? 2 - p 2 - m 2 N + it v 



7°-p7 + m N ) 



Since t — > +00 then t — x' a > and we close the in- 
tegration contour over R° with a semi-circle of infinite 
radius on the lower half-plane picking up the pole at 
R = p° — it = E(jp) — it. Applying the Dirac equation 
to the result one arrives to eq.(||). 

Then taking into account eq.mh and the expansion 
D-V = [D - A]- 1 ^ = DjV + D^AD^-f + .... 
we can rewrite eq.(0) as: 



N 

: J [dU](det D) e ^ dxC ™ £ ( CT ) II / dlin J dx ^ 



'^(Pn) 



8(Xn - X CT „)(5„ 



i dt /df 



A{x n ,x an )e- iE ^ t ' -i dt dt' dzdz'e^^ 



A(x n ,z)D- 1 (z, z')A(z', XfT Je- iE ^)t' + 
e^^u(p a J 



(6) 



where now both t and t' are integration variables being 
the time components of x n and x Gn , respectively. The 
dots just refer to those terms with an increasing number 
of insertions of the operator A coming from the geo- 
metric expansion of D^ 1 = [D — A] 1 discussed above. 
Furthermore A(x, y) is defined such that 



dxip{x)A{x)ip{x) = J dxdyif)(x)A(x,y)ip(y). 
As a result of cq.(|^) we can simply state that: 

[dU] (det D) e i $ dxC ™ (det T) (7) 



iv.a.s,', 



where the tilde in det T indicates that the determinant 
has to be taken in the subspace of the Fermi sea states 
expanded by the basis functions e lp " x u(p„) with |p n | < 
hp. In this notation T is given by: 



T = h- il dt dt 



JdtJ dt' e lHot ~/°A [h - D^A] 1 e " 



iH t' 



where I3 = 5(x„ — x m )<5„ im and analogously I4 = 6(x r , 
x m )8 n ,m- On the other hand, e~ lHot e jp " x u(p n ) 



-ip n x 



'u(pn) and e lp " x u t (p„)e 



iH t _ „ip n x.A 



W T (P")- 



In order to obtain from eq.(|7j) the contributions of 
the surrounding medium to the generating functional 
it is convenient to exponentiate det T as exp(Tr log.? 7 ) 
(where the tilde has the same meaning as before). Then 
we have: 



e iz[v, a ,s, P ] = j [ d[/ ]( det £>) exp L I dx £^ + 

•* ^ ^ r=l a=l 



e ipy u ra (p)\ (8) 



where we have indicated explicitly the spin r and isospin 
a indices. It is very appropriate to stress at this point 
that eq.(^), although formal, is non-perturbative. From 
this result we can readily read out the new effective chi- 
ral Lagrangian density C wv in the presence of a nuclear 
density just by equating the expression between curly 
brackets to ijdx C„ T . 

3. The perturbative theory is obtained by expanding 
log J 7 in eq.(||), so that e tZ can be written as: 

J[dU] (det D) exp j ijdx £„„ - % £ E / {2v)nE(p) 



dxdye tpx u ra {p)A[I A -D l A] \ xa . ya) e lpy u ra (p) 



1 j2 f kF dp f kF dq 



dx dx' dy dy 



2 r% a 7iiJ (27r) 3 2E(p)7 {2^f2E{q). 
e ipx u ra (p) A[h - D~ 1 A}~ 1 \ [xatXlf3) e-^'usp (q) e^' u sf3 (q) 

A[h - D^A]-\ y ,^ y<x) e-™u ra (p) + ...j. 

2 

Finally, taking into account the relation ^ u r (p) (8> 

r=l 

u r (p) =]/ +rnjv we arrive at the following expression 
for the generating functional: 



e lZ = J [dU] (det D) exp i^ij dx - i J- 



kF dp 



(2tt) 3 2E( P ) 

dx dy e lp[x - y) Tr I A[h - D^A]' 1 \ {x<y) +m N ) ) + 



1 f kp dp 



dq 



dx dx' dv dv' e ip{ - x ~ v ^ 
2 J {2nf2E{p)J (2nf2E{q)J * yay 6 

e-^'-v') TrU[J 4 - D^A]-\ XiXl) ^+m N )A[h 
DQ X A]-\ v , ty) (tf+m N )\ + ... \ 



(9) 



where the trace refers both to the isospin and spinor 
indices. The previous formula implies a double ex- 
pansion for obtaining the contributions of the nuclear 
medium to the chiral Lagrangian. One is the stan- 
dard chiral expansion by expanding the vacuum opera- 
tor A[h - D^A}- 1 = A + AD q 1 A + AD q 1 ADq 1 A + ... 
together with A itself, A = A^ + A^ + A^ + ... in 



3 



increasing powers of momenta and quark masses valid 
at low energies. Explicit expressions of A up to 0(p 3 ) 
can be obtained from ref. The another is an ex- 
pansion in the number of insertions of on-shell fermions 
belonging to the Fermi sea, schematically indicated in 
fig.l by a thick solid line. Both expansions can be re- 
lated by giving a chiral power couting to kp which can 
be naturally counted as 0{p) [M since for nuclear satu- 
ration density kp ~ 2M W with the pion mass. More- 
over, the circles labeled by T correspond to the non-local 
operator —iA[I± — Dq 1 A]~ 1 . Note that when inserting 
n Fermi seas from the expansion of the logarithm one 
picks up a factor — (— l) n /n where the global minus sign 
appears due to the fermionic closed loop, n is a com- 
binatoric factor because any cyclic permutation in the 
trace of n Fermi seas with their associated n T operators 
gives the same result and finally the sign (—1)™ is a pure 
in-medium factor that one has to keep in mind and is 
already present in standard many-body theory JtJ . 



r 1 





r ) ( r ) 




FIG. 1. Diagrammatic expansion of eq.(|9|). Every thick 
line corresponds to the insertion of a Fermi-sea and each circle 
to the insertion of an operator Y = —iA [I4 — D^A] 

A generally non-local vertex T comes from the itera- 
tion of the A operator with intermmediate free baryon 
propagators Dq 1 , obeying the usual Feynman rules, see 
fig. 2. Notice that A was defined from the Lagrangain 
rpDrp removing the free term tpD^ and changing the 
sign to the rest, this is why a minus sign appears in 
front of A in fig. 2. On the other hand i Dq (x,y) — 



i f d p 



is the usual baryon propagator. 



FIG. 2. Expansion of the generalized non-local vacuum 
vertex F. Every solid line corresponds to a vacuum baryon 
propagator and each circle to the insertion of an operator 
— iA from ^pDip. 

Hence a final diagram, when expanding T up to the 
required accuracy, will be a set of n > 1 Fermi-sea in- 
sertions, m > free baryon-propagators and of m + n 
vertices — iA. First we include the —1 global sign because 
of the fermionic closed loop, and the combinatoric factor 
1/n together with the sign (—1)". Then, following the 
diagram in the opposite sense to that of the fermionic ar- 
rows, write for each Fermi-sea an integral J kF 
with p° = E(p), for each vaccum baryon propagator 
with free momentum p write i f 2 ^ +m i N , - and for 

a vertex in momentum space a term — iA(27r) 4 , keep- 
ing in mind the energy-momentum conservation at each 



vertex. Finally sum over the spin and isospin indices of 
the fermions. This defines explicitely the Feyman rules 
in momentum space in order to obtain i(2-7r) 4 times the 
desired connected graph accompied by the global delta 
of energy-momentum conservation. Analogous Feyn- 
man rules hold of course in configuration space, e.g., see 
eq-(|). 

An equivalent way to state the previous rules, without 
including explicitely the integral symbols and factors 27r, 
is to write the same sign-combinatoric factor — (— l)"/n 
and vertices — iA as before. Then for the free nucle- 



< 2 —m 2 +«e 



and for 



ons one has just the free propagator i- 
the Fermi-sea baryons the factor {2ir)8(p 2 — m 2 )8(p°)(p' 
+m)8(kp — |p|). Finally sum over all the discrete indices 
attached to the fermions and integrate over all the free 
four-momenta with the measure J d 4 p/ (2n) 4 after tak- 
ing into account energy-momentum conservation at each 
vertex. 

Fig.l fixes the skeleton structure of standard in- 
medium vertices since still one has to consider the pion 
fields contained in A over which one has to integrate in 
eqs.(||), (H) to finally obtain the generating functional. 
That is, from the vertices A as well as from C^, one 
can generate internal as well as external (coupled to the 
sources) pionic legs denoted by dashed lines in figs. 3a and 
3b. Here one has essentially the same Feynamn rules 
than in vaccum in order to proceed in a perturbative 
way, simply for each pionic line with four-momentum q, 

for 



(2;r) 4 q' 2 -M' 2 +ie ' 



one writes the vacuum propagator i J 
the first version of the in-medium Feynman rules. For 
the second one has 



q ^ -M 2 +it • ^ ne important remark to 
keep in mind is that a generalized vertex has analogous 
properties to those of a standard local quantum field the- 
ory one to the effects of determing the numerical factors 
accompanying the exchange of pion lines inside a given 
diagram. This can be seen just by applying standard 
perturbative techniques in path integrals to the action 
given between brackets in eq.(^J) or more explicitely in 
eq.(||). Several examples are discussed in detail in ref. 

B 

4. Notice also that ultra-violet parts in the integra- 
tion over running pionic momenta generate local multi- 
nucleon vertices and hence a proper treatment of pion 
loops can only be done by including simultaneously local 
nucleon interactions in order to reabsorb the divergences. 
This is denoted by the squares in figs. 3c and 3d. It is 
important to realize that while in fig. 3a the pion is ex- 
changed inside the same vertex (in the generalized sense 
of above) in fig. 3b the pions are exchanged between two 
of them. This implies the presence of only one Dirac 
trace in fig. 3c (the flow of the Dirac indices along the 
propagators on both sides of the square is indicated by 
the arrows of the open solid lines) and two Dirac traces 
in fig. 3d. Thus, when writing the matrix elements corre- 
sponding to diagrams of the type of figs. 3c and 3d (with 



4 



local interactions) one has to keep track of the number of 
closed traces in the spinor indices since each of them will 
lead to a sign (—1) and its own combinatoric factor 1/rij 
(with rij the number of Fermi-sea insertions in the closed 
Dirac loop) together with the sign (— l) nj . One can also 
arrive to the same conclusions by considering that the 
dashed lines originate from the exchange of an arbitrary 
heavy meson (the only relevant point in our derivations 
is the bilinear character of the meson- N N interaction 
in the nucleon fields). In this way it is straightforward 
to realize the presence of a factor 1/2 in front of fig. 3d 
together with the factors ^) nj+1 / n j which here is 

simply 1. The relation between contact four nucleon in- 
teractions at low energies in effective field theories and 
its saturation by integrating-out heavy meson resonances 
has been established in ref. Jlal. 



momentum regime. However, when a Fermi-sea baryon 
line is present the momentum running though the baryon 
lines is of the form pj = p + Qj with p = (E(p), p) and 
|p| < hp. Thus one is inserting a medium parameter Uf 
and shifting upward up to around tun the energy level 
around which one is perturbating. 

5. We now turn to the generalization of the formalism 
to the case of asymmetric nuclear matter, with different 
densities of neutrons, p n , and protons, p p , with Fermi 

momenta k ( p ] = (37r 2 p„) 1/3 and fc$ = (i^pp) 1 ^ 3 , re- 
spectively. Following the previous derivation of eq.(||) 
one can easily convince oneself that the only change is 
to remove the sum over isospin indices and to distinguish 
between a = 1 (proton) and a = 2 (neutron). In this 
way we will have: 



a) c) e) /■** dp t f . t , N , 





2 rk[ 



b) 



d) 



FIG. 3. Diagrams 3a and 3b represent some typical 
many-particle diagrams generated from eqs.Q, The cir- 
cle indicates an A operator insertion (which in addition can 
have attached to it more lines than shown) and the dashed 
line corresponds to a pion exchange. Figs. 3c and 3d arise by 
considering the ultra- vioalet divergent part of the pion loops 
leading to local terms denoted by squares. Fig.3e is a local 
7riV counterterm, indicated by a diamond. For more details 
see the text. 

There is still an important difference to be discussed 
when comparing figs. 3c and 3d, which in fact is related 
to the presence of the factor (det D) in all the formulae 
from cq.(^) to eq. ([)]). The latter corresponds to contri- 
butions to the chiral Lagragian from closed fermion loops 
in the vacuum and in the spirit of the effective field theo- 
ries these contributions, coming from states with masses 
close to or above the chiral scale A x rj M p , are incorpo- 
rated in the counterterms of the vacuum effectrive field 
theory. In this way we will set det D = 1 in the fol- 
lowing. Indeed on the right hand side of the square of 
fig. 3c we can recognize a momentum loop flowing along 
the free baryon propagators without any insertion of a 
Fermi-sea. This simply means that since all the mo- 
menta that could go into this closed momentum loop 
are soft because they come from the circles, it just cor- 
responds to vacuum renormalization from heavy parti- 
cles and is reabsorbed in the low energy counterterms. 
This is schematically shown in fig.3e with a diamond cor- 
responding to a higher order ttN counterterm. This is 
consistent as far as one is restricted to the low energy and 



uU(p) log[^]| (x2 , y2) e™u r2 (p)\ 



(10) 



Notice that we have indicated separately the ener- 
gies of protons E(p)i and neutrons E(p)2 with three- 
momentum p, since the previous equation is valid for 
the non-isospin limit as well. Nevertheless, in order to 
simplify the formulae, we will consider in the following 
the case with equal nucleon masses. We also introduce 
the 2x2 matrix: 



n(p) = 









0(k { F ] -p) 



( nip)! 

n(p) 2 
n(p)I 2 + n(p)r 3 



(11) 



with I2 the 2x2 unity matrix, T3 the usual Pauli ma- 
trix diag(l,— 1), n(p) — (n(p)i + n(p)2)/2 and h(p) = 



(Pk- 

[n{p)\ — n(p)2)/2. Then eq.([10|) can be rewritten as: 
e«^.,] = J m exp {, j^j^jdyj^-^ 
e -P(x-y) Tr ^ n{p) log[jr] 1 (x y) ( y +m) 

dp 



[dU]e^^JdxC^-ij {27r)32E{p) 
Tr[A[I 4 - D^A\-\ X ^ +m N ) n(p) ) + 



dx dy e 



ip(x-y) 



5 



1 f dp 



r/q 



dx dx' dv dv' e lp< - x ~~ v ^ 
e-^'-y'^Tr^Aih-D^Aj-^^W+m^niq) 
A[h-DQ X A]-\ y , >y) {tf+m N )n{p) ] + ... L. 



(12) 



Comparing this equation with eq.(^), the only difference 
in the Feynman rules is the inclusion of an isospin matrix 
n(p) associated to every Fermi-sea insertion. 

The case of a polarized nuclear matter can be treated 
in the same way just by doing in eq.(|i"c|) the replacement: 

/ d P U t a (P) -■ "™(p) -» / rfp4ia(p) -■ u r ia (p) 
dP*4 2 a(P) •■• U r2a {p) (13) 

since there are two spin states per momentum state. 

6. The present formalism is applied in rcf. |l4|] to eval- 
uating several quantities relevant for low energy QCD 
in the nuclear medium. There we will also address in 
detail the issue of the chiral counting from eq.fll^) in 
the nuclear medium and the limitations of a plain per- 
turbative treatment of CHPT at finite density. Just to 
mention that instead of a relativistic treatment of the 
baryons we could also have considered in the same way 
the non-relativistic case. This makes more straightfor- 
ward the evaluation of baryon loops [^6) although one 
has also to take into account recent developments in the 
field of effective field theories with propagation of rel- 
ativistic heavy particles |l7 19 1. In any case the main 
issue still to be addressed in the medium, as discussed 
in refs. 16 1^], is to deal with the problem of the large 
S-wave scattering lengths in the nucleon-nucleon scatter- 
ing which introduces a new extra scale of only ~10 MeV 
alredy in the vaccum case, where consistent power count- 
ing schemes have been developed 20 . Some interesting 



findings in this direction, requiring further consideration, 
can be found in rcf. [pl| for a theory without pions. 

To conclude, we have derived the SU(2) chiral La- 
grangian with extenal sources in the presence of non- 
zero nuclear density by explicitely working out in quan- 
tum field theory the in-medium contributions from the 
ttN chiral Lagrangian, eq. (^) (symmetric nuclear mat- 
ter), eq.(|l^) (asymmetric nuclear matter) and eq.([l3|) 
(asymmetric and unpolarized nuclear matter). Then the 
perturbative theory is developed and the corresponding 
Feynman rules are given. Contrarily to standard many- 
body techniques, the rules and diagrams here derived for 
the general relativistc case are analogous to the usual 
ones from vacuum quantum field theory without modifi- 
cation of the baryon propagators establishing a neat sep- 
aration between in-medium and vacuum contributions, 



as indicated in figs.l and 2. Applications of this many- 
particle relativistic quantum field theory formalism to 
actual calculations can be found in ref. O]. 
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